2, 2008 6:7 WSPC/INSTRUCTION FILE mgb2'MPLB 



POINT-CONTACT SPECTROSCOPY OF MAGNESIUM DIBORIDE 
WITH DIFFERENT COUNTER-ELECTRODES 



A. I. D'YACHENKOt, V. YU. TARENKOVt, M. A. BELOGOLOVSKIlt, 
V. N. VARYUKHINt, A. V. ABAL'OSHEV** AND S. J. LEWANDOWSKI* 

t Donetsk Physical and Technical Institute, National Academy of Sciences of Ukraine, R. 
Luxemburg 72, 83114 Donetsk, UKRAINE 

^Instytut Fizyki Polskiej Akademii Nauk, Al. Lotnikow 32, 02-668 Warszawa, POLAND 

abala&ifpan. edu.pl 

Received February 2, 2008 

We report on tunneling and Andreev-reflection conductance spectra of 39 K supercon- 
ducting magnesium diboride, obtained with Pb and Au counter-electrodes. Two distinct 
steps at close to 2.7 and 7.1 meV appear in a low-resistance metallic-type Au-MgB2 
junction characteristic, whereas a tunneling-like spectrum measured for the same junc- 
tion, annealed by the application of dc current, exhibits only a rounded contribution 
of the larger gap. Junctions with a superconducting lead counter-electrode pressed into 
a bulk MgB2 sample reveal two conductance peaks that are interpreted as the result 
of the formation of highly-transmitting break junctions inside the magnesium diboride 
ceramic. Our results strongly support the two-band model with two different gap values 
on quasi-two-dimensional cr (7.1 meV) and three-dimensional n (2.7 meV) Fermi surface 
sheets of MgB2. 

PACS Number(s): 74.25.Kc, 74.45.-l-c, 74.50.+r, 74.70.Ad 
1. Introduction 

The discovery of superconductivity in the simple binary compound, MgB2, gen- 
erated an intense research activity - mainly, because MgB2 exhibits a transition 
temperature Tc of 39 K, which is almost two times higher than that of any other 
intermetallic superconductor known earlier ^ , but also because of good perspectives 
for practical applications. During the two years, which have passed from this dis- 
covery, our understanding of physical properties of the novel superconductor has 
made rapid progress. 

After confirmation ^ of the spin-singlet type of the Cooper pairs involved, the 
main discussions relating to fundamental physics concentrated on two issues: the 
mechanism of the Cooper pair formation and the symmetry of the order param- 
eter. The boron isotope effect experiments ^ were the first indication of phonon- 
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mediated superconductivity. Direct experimental probes of phonon-induced struc- 
tures in point-contact ^'^'^■'^ and tunneling ^ superconducting junctions have pro- 
vided unequivocal support to the theories f .10,11,12,13,14^ which explain the relatively 
high Tc as originating from the strong anharmonic electron-phonon coupling to 
Raman-active high-frequency E2g vibration modes, corresponding to the in-plane 
distortions of the boron hexagons. 

It took much longer to elucidate the properties of the superconducting order 
parameter, but now from the results of tunneling and photoemission spectroscopy, 
specific heat, and other measurements emerged strong experimental evidence that 
magnesium diboride belongs to a rare class of materials with two gaps of different 
widths, which arise from a huge anisotropy of electronic properties (c/the review 
of tunneling data for magnesium diboride ^^). The Fermi surface of MgB2 is very 
anisotropic and consists of four bands: two cr-type quasi-two-dimensional cylindrical 
hole sheets and two 7r-typc three-dimensional tubular networks, exhibiting very 
different electron-phonon coupling strengths As a result, the superconducting 
gap varies considerably on the Fermi surface and, as follows from first-principles' 
calculations ^'^^ clusters into two groups. The a cylindrical sheets relate to a larger 
gap of an average value of 7 meV with small variations, while the smaller gap ranges 
from 1.2 to 3.7 meV on the three-dimensional tt sheets. These calculations support 
a two-gap scenario for magnesium diboride, originally proposed by Shulga et al. 
in order to explain the behavior of the upper critical magnetic field. It should 
be noted that, in contrast to conventional superconductors, impurity scattering 
does not average out strongly different gap values. Averaging becomes ineffective 
in MgB2 because the different symmetry of a- and 7r-type bands makes interband 
scattering much weaker than intraband impurity mixing of electron states 

Based on such ideas, an effective two-band model with ^„{T = 0) = 7.1 meV 
and A^(T = 0) = 2.7 meV was derived, and consequences for tunneling into MgB2 
were deduced The discussion of our tunneling and point contact data relating 
to the energy gap distribution will be mainly based on this model and will use a 
minimal number of fitting parameters. Most samples studied previously were in the 
intermediate regime between a direct metallic contact and a tunnc;! junction with a 
strong insulating (I) barrier. The main conclusions for these structures have been 
drawn by fitting the differential conductance data to the one-dimensional Blonder- 
Tinkham-Klapwijk (BTK) model modified by introducing a damping parameter 
r, and by taking into account two possible channels for the charge transport, so that 
not only the gap values, but also the ratio of the individual channel contributions 
were taken as adjustable parameters. In this paper, we present the data for low- 
resistance junctions with a gold counter-electrode, which to our knowledge are the 
first reported direct metallic- type contacts with superconducting MgB2. and which 
can be described within the two-band model ^'^ without any additional parameters. 
Moreover, by applying dc current pulses we have successfully transformed the clean 
contact into a high-resistance tunneling junction, and reproduced in this way the 
two limiting regimes of the theory The second issue, which is discussed in the 
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paper, is the assertion by Schmidt et al. ^^'^^ that pressing a counter-electrode into 
the MgB2 sample often breaks off an MgB2 crystal fragment and results (even if 
wc arc dealing with a normal (N) injector) in a junction between two (identical) 
superconductors, instead of the expected normal metal - superconductor structure. 
Such internal break-junctions apparently can be formed also when a superconduct- 
ing counter-electrode is pressed into MgB2. Our experiments with a lead injector 
provide an independent confirmation of the observation made by Schmidt et al., 
and fit into the double-gap scenario as well. 



2. Results and Discussion 
2.1. Sample preparation 

We used polycrystalline MgB2 samples in the form of 15 x 1 x 0.1 mm^ rectangu- 
lar bars of high density, prepared by compacting commercial magnesium diboride 
powder at ~ 20 — 30 kbar. Magnetization measurements have shown a supercon- 
ducting transition with the onset and midpoint at 39.0 and 36.5 K, respectively. 
All experiments were performed at 4.2 K. The junctions based on these samples 
exhibited mostly no 'native' barrier, and their resistance varied from several to a 
few tens of Ohms. In order to establish a suitable barrier between the electrodes, 
we applied short current pulses with an amplitude of several tens of milliamperes 
and a duration of a few seconds; as a result we obtained junction resistances up 
to a few hundred Ohms. Two different counter-electrodes have been used: lead and 
gold. Current / versus voltage V characteristics as well as differential conductance 
spectra G{V) = dI{V)/dV were measured by conventional four-probe lock-in tech- 
nique In the paper we present typical data obtained in this manner. 

The junctions with Pb injector were fabricated by mounting a Pb hemisphere 
of 0.5 mm diameter on an elastic beryllium strip, which was driven into the MgB2 
surface. Point-contact junctions with an Au counter-electrode were prepared by 
touching the superconducting surface with a 0.02 mm gold wire wrapped around 
the MgB2 sample; the formation of SNS junctions was avoided in this manner. The 
device was stabilized by a coat of varnish. This procedure, which yielded junction 
resistance of about 10 17 and conductance spectra indicating a direct heterostruc- 
ture formed by normal and superconducting electrodes, could in principle produce 
several point contacts instead of a single one. Such occurrence is not important 
in the case of clean junctions, since in such a case the spectral contribution of 
the contacts with transparency significantly smaller than unity is negligible. In our 
best Au-MgB2 samples, an almost ideal metallic Sharvin contact between gold and 
magnesium diboride was formed. This can be seen from the conductance drop with 
voltage increase from y = to A^/e (figure 2). The drop is almost by a factor of 
two, with the small deviation attributed to damping and thermal smearing effects. 
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2.2. Pb-MgB-i junctions 

The junctions with a Pb injector were prepared in an analogous manner to that used 
by Schmidt et al. , who formed symmetric MgB2 — MgB2 structures by breaking 

off an MgB2 crystal fragment under the tip pressure. Following the arguments of 
the cited work, we also suppose that our Pb-MgB2 samples are symmetric break 
junctions inside the magnesium diboride bulk, rather than contacts between two 
different superconductors. The nature of the junction barrier is unknown, but its 
small resistance compared to that of the SIS structures suggests that we are deal- 
ing rather with SNS contacts (see also the paper of Gonnelli et al. and arguments 
given below). The initial liV) characteristic exhibited a critical supercurrent at zero 
bias /c ^ 0.9 mA and normal-state resistance i?N = 1-8 ri, i.e., Jci?N ^ 1-6 mV [plot 
1 in the inset in figure 1(a)]. This value of /ci?N is significantly lower than that pre- 
dicted by the model ^® for a Josephson current across a symmetric MgB2 — MgB2 
device, but is consistent with /c-Rn values found previously by two other groups ^^'^^ 
for break junctions of magnesium diboride. 

The evolution of the interface properties of the junctions, caused by the appli- 
cation of the forming current pulses, is clearly visible in the current-voltage char- 
acteristics shown by the inset in figure 1(a). After two successive current pulses 
(characteristics '2' and '3' in the inset), we obtained a finite resistance i?N = 3.6 f2 
at zero bias (characteristic '3'). The latter characteristic, shown in the main panel 
of figure 1(a) in a larger voltage scale, exhibits an excess current /exc known to 
be the manifestation of Andrccv reflections in high-transparency superconducting 
heterostructures. Remarkably, only twofold increase of junction resistance was suf- 
ficient to suppress the supercurrent, whereas the Josephson current in reference 
was visible up to contact resistances of about 10 kf7. Our result can indicate that we 
create multiple incoherent-scattering centers and the Josephson current is cancelled 
because of the phase-breaking events. Liao et al. found that magnesium diboride 
samples prepared by solid-state reaction contain large amounts of precipitates com- 
posed mainly of magnesium and oxygen. We believe that their presence is the origin 
of a strong damping factor in conductance spectra of MgB2 (even without any cur- 
rent treatment), stressed in almost all publications on point-contact spectroscopy 
of this compound. It is highly probable that by applying currcint pulses, we create 
such centers and destroy the phase coherence during the charge transfer across the 
interface region between two superconducting banks. 

To further verify the hypothesis of the SNS symmetric break-junctions formed 
inside MgB2 bulk, let us analyze the high-voltage asymptote of the current-voltage 
characteristic shown in figure 1(a) and known to be determined by a set of inter- 
layer transparencies Their distribution has been found in two extreme cases: 
for a diffusive normal conductor by Dorokhov and for a disordered insulating 
interface by Schep and Bauer (see also Naveh et al. ^^). For two identical super- 
conductors both theories predict I{V) = G-^iV + fcA/e), where the last term is the 
excess current /exc- For the Dorokhov distribution, the coefficient k = 1.467, while 
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Fig. 1. Current and differential conductance in function of voltage for a junction formed by press- 
ing Pb into MgB2 bulk sample: (a) I — V characteristic at T = 4.2 K of a current modified 
junction (marked 3 in the inset) exhibiting an excess current /exc = 1-6 mA (solid line in the 
main panel); the dotted line was measured in the normal state. The inset shows successive mod- 
ifications of the I — V characteristics (labelled 1 to 3) after applying short current pulses, (b) 
differential conductance normalized to the conductance at 20 mV corresponding to characteristic 
3 demonstrates a prominent zero-bias maximum with two peaks on each side (arrows indicate 
voltage positions of the combinations of small and large gap values 2A^/e and Aa-/e + A^/e). 
The inset shows G(V')/Gn corresponding to characteristic 2 as a function of A/eV (continu- 
ous line) and the plot of the theoretical -^^ interpolating formula for disordered SNS junctions 
G(V)/Gn = 1 + 0.5i/A/ey (dashed line); the best fit was found for A = 6.8 meV. 

within the Schep-Bauer approximation, k = 1.055 '^'^ . Applying the arguments of 
reference , we presume that for dissimilar superconductors with gap values A and 
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A' the excess current in a SNS' junction is given by /exc^N ~ A(A + A')/(2e). We 
can perforin then an analysis similar to that presented in reference ^® for the dc 
Josephson current. For a symmetric MgB2 — MgB2 structure we conclude that the 
Jexc-RN product should vary from < 4 mV (tunneling along the c-axis direction) to 
6 mV (tunneling within the a — b plane) in the case of a S-diffusive N-S structure, 
and between 3 mV to 5 mV in the case of a disordered insulating interface. At the 
same time, for a contact formed by a conventional superconductor like Pb with a 
single gap equal to 1.4 mV, the /exc-RN value cannot exceed 4mV in any model. 
From figure 1(a) we find /exc-^N = 5.7 mV, i.e. a value much greater than the latter 
limit. This strongly supports the thesis that the weak link in our structures with 
lead counter-electrode is not at the interface of the contacting materials but rather 
inside the magnesium diboride bulk. Furthermore, the product /exc^^N exceeds the 
upper limit for the Schep-Bauer distribution, providing an argument for the SNS 
nature of the junction. 

To verify further this scenario, let us compare our dc differential conductance 
with the predictions of the scattering-matrix theory of multiple Andreev reflections 
in short disordered SNS junctions Using zero voltage asymptotic behavior of the 
dc differential conductance G{y) = dIiy)/dV (equation (14) of reference ^'^), we 
can check on the background behavior of our conductance spectra. We apply an 
interpolating formula G{V) = GN(l + 0.5tanh(A/2T) y^A/eV), which is vahd both 
at very low bias (the near-zero-bias 1 / ^/V singularity is the hallmark of Andreev- 
reflection processes occurring at each junction face), and at very high voltages (in 
the latter case GiV) is equal to the normal-state value Gn = I/^^n)- Using A as 
a fitting parameter, we found that the best agreement with the experimental data 
is achieved for A = 6.8 meV. The overall behavior of the experimental normalized 
conductance G(T^)/G'n corresponding to the experimental data and plotted in func- 
tion of ^6.8/ey is satisfactorily described [see the inset in figure 1(a)] by a simple 
fitting line G(F)/Gn = 1 + CS-y/A/eV in the voltage range above 0.5 mV, where 
temperature smearing is not effective. In the small bias region up to 15 meV, the 
conductance exhibits a fine structure, which according to Bardas and Averin can 
be attributed to subharmonic gap singularities, as well as to the gap features. 

Characteristic '3' is more smooth than '2', and the corresponding differential 
conductance presented in the main panel of figure 1(b) contains a strongly sup- 
pressed zero-bias peak and two maxima. Their voltage positions are very close to 
the combinations of small and large gap values 2Ao-/e and Ac^/e + A^/e expected 
for the magnesium diboride and shown in figure 1(b) by arrows. We should em- 
phasize that we do not observe any signs of the lead gap. Moreover, the features 
discussed above, including the peak at F = 0, are seen at temperatures exceeding 
the Tc of Pb, and this rules out their possible relation to the lead counter-electrode. 
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2.3. Au-MgB-i point contact junctions 

Let us examine now the conductance spectra of our Au-MgB2 point contacts in 
the framework of a double-gap model with fixed gap values and weighing fac- 
tors These quantities arc found from first-principles' calculations, therefore the 
adjustable parameters are reduced to the barrier strength and the damping factor 
r. The latter one was introduced into the tunneling density of states by Dynes et 
al. and is often referred to as the Dynes factor. The barrier strength is mea- 

d 

sured by a parameter Z = J H{x) dx/hvp where H{x) is the barrier height, 



X denotes the coordinate normal to the interface, d is the barrier width, and vp 
is the x-component of the Fermi velocity. The conductance spectrum G{V) is as- 
sumed to be a sum of two single-mode BTK contributions related to the two 
a- and 7r-type bands in proportions strongly depending on the injection angle. The 
relevant weighting factors, which determine the ratio of the two contributions, are 
proportional to the squared plasma frequencies and are given in reference for the 
transport within the a — 6 plane, Ga-b{V), and for tunneling in the c-axis direction, 
Gc{V). In the present experiment, the conductance is a mixture of a — 6 plane and 
c-axis contributions with unknown coefficients, which can be different for different 
polycrystalline samples. Only to get a qualitative estimate of the averaged conduc- 
tance spectra predicted by the theory, we compare the measured spectra with the 
mean characteristic Gm{V) = [Ga-b{V) + Gc{V)]/2. As it was stated above, we 
prefer to obtain qualitative conclusions with a limited number of fitting parameters 
rather than to introduce weighting factors with the aim to improve the agreement 
between theory and experiment. 

In reference , a two-step structure was predicted for a direct contact between 
a normal-metal injector and magnesium diboride and for charge transport in the 
a — b plane (upper curve in the inset in figure 2 calculated for T = 0). For the c-axis 
transport, the main contribution comes from the tt band (curve 2 in the inset). The 
averaged characteristic (see curve 3 in the inset and the dashed curve in the main 
panel), which takes into account a small Dynes factor F and temperature smearing 
effects determined by a standard procedure involving Fermi distribution functions, 
qualitatively well reproduces the two-step feature observed by us expca'imcntally 
(solid line in figure 2). It should be stressed that this good agreement of experiment 
and theory was obtained with a single fitting parameter F. 

We formed next a low-transparency insulating (I) barrier in our Au-I-MgB2 
junctions by applying additional current pulses and measured the resultant G{V) 
characteristics in a wide range of bias voltages (up to 100 meV). The conductance 
background usually exhibited a roughly parabolic dependence on voltage, with a 
very small offset of minimum from V = (sec figure 3), indicating an almost sym- 
metrical barrier and proving the tunneling character of the charge transport. A 
simple formula for a rectangular potential barrier may be used to find an average 
barrier height H and its thickness d. In order to carry out the evaluation, we need 
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Fig. 2. Conductance spectrum at T = 4.2 K of a direct Au-MgB2 contact (i?N = 20.5 f2), normal- 
ized to the conductance at 20 mV (solid line) . The dashed line represents a mean characteristic 
G.m(y) calculated from the two-gap model in the limit of a completely transparent interface 
and Z = 0, small smearing parameter T = 0.23 and T = 4.2 K. The inset shows zero-temperature 
theoretical plots of conductance {Z = 0, T = 0): 1 - charge transport within the a&-plane, 2 - 
transport along the c-axis, 3 - mean conductance Gm{V); the plots are displaced vertically for 
clarity. 




Fig. 3. Conductance at T = 4.2 K of a tunneling Au-l-MgB2 junction with normal-state resistance 
Rn = 1.2 kfi formed by applying several current pulses (solid line). The dashed line shows a 
parabolic fit. 
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the contact area, which is unknown. Nevertheless, a minimal value for Z parameter 
can be estimated from the parabolic shape of the conductance curve. According 

to reference "^^j the coefficient k in the quadratic term in the voltage dependence 
of conductance, G{V)/G{0) = 1 + [eV/nY, is given by k = H^'"^ / {Q.lS2d), where 
H is in volts and d is in Angstroms. We assume d to be not less than two- three 
lattice parameters and take as a minimum d = 10 A. Substituting n found from the 
measured conductance background (figure 3), we obtain Z > 0.85. Consequently, 
we take the first fitting parameter Z equal to 0.9. 




Voltage (mV) 

Fig. 4. Tunneling-like conductance spectrum of Au-I-MgB2 junction with normal-state resistance 
Rn = 1-2 kQ normalized to the conductance at 20 mV (solid line). The dashed line represents 

the calculated mean characteristic Gm{V) following the two-gap model for Z = 0.9, F = 2.9, 
T = 4.2 K. The inset shows zero-temperature theoretical curves for Z = 0.9 and vanishing damping 
parameter F: 1 — charge transport within the afe-plane, 2 — transport along the c-axis, 3 — mean 
characteristic Gm(V). The plots in the inset are displaced vertically for clarity. 

In figure 4, the main panel presents the measured conductance spectrum of the 
same junction in the near-gap voltage region, while in the inset we show the con- 
ductance calculated for tunneling into the a—b plane and along the c-axis direction 

for vanishing Dynes factor, F = 0. The latter characteristics clearly indicate that 
the most pronounced feature in the average conductance spectrum should be the 
smaller gap peak, whose signature is missing from the measured spectrum. The 
most probable cause is the strong disorder induced by the forming current pulses 
and the resulting significant inelastic scattering near the interface. For large values 
of F (the second fitting parameter) the agreement between the experimental and 
calculated conductances is very good, as shown in the main panel in figure 4. 

To conclude the discussion of near-gap features in the junction characteristics, 
let us make two remarks concerning the role of disordering effects. The analysis 
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based on a simple BTK approach for describing the experiment with (probably) 
numerous elastic and inelastic scatterers randomly distributed within the barrier 
appears to be somewhat naive, nevertheless is jiistified, since one can always in- 
troduce an effective transparency parameter Z within a single-mode model by 
integrating over a certain distribution of transparencies (see the corresponding dis- 
cussion in reference ^^). It should be stressed that the additional averaging already 
takes into account the possible appearance of multiple point contacts in our experi- 
ments. The averaging provides some additional smearing, but the most transparent 
contacts will dominate in the spectrum. 

The second remark concerns the possibility of observing double-gap spectrum in 
a disordered superconductor with important charge scattering. Our samples contain 
comparatively large grains. Therefore, we expect that the measured characteristics 
of a clean point contact reflect the bulk properties of the material. It is known that 
dense bulk samples of MgB2 are usually in the clean limit, with the elastic mean free 
path large in comparison to the superconducting coherence length. This condition 
is necessary to observe two distinct gaps in a sufficiently clean superconductor. At 
the same time, our method of fabricating the insulating layer between the electrodes 
generates, as it was argued above, a great amount of inelastic-scattering inclusions. 
At first sight, the interband impurity scattering should cause the gaps to converge. 
But this is not the case for magnesium diboride. According to reference , optical 
properties of dense polycrystalline samples of MgB2, as well as a surprisingly small 
correlation between the defect concentration and T^, can be understood only by 
assuming negligibly small inter-band impurity scattering. The possible reason is in 
the specific electronic structure of MgB2 , namely the the two-dimensional cr-band 
states overlap slightly with the states in the Mg plane, where defects are most likely 
to occur ^® . Thus one can expect that the interband scattering rate between the a- 
and TT-bands is small, even in low-quality samples. 



3. Conclusions 

In summary, we have realized several different types of heterostructures based on su- 
perconducting magnesium diboride: from a low-resistance (almost ideal) Au-MgB2 
point contact to an SNS junction (with a lead counter-electrode) obtained by apply- 
ing current pulses. Moreover, with this technique we have succeeded to transform 
point-contact samples with gold injector into tunnel-like devices. In the case of lead 
electrode, we have confirmed the observation of Schmidt et al. that pressing the 
electrode can result in the formation of symmetric superconductor-superconductor 
junctions inside the magnesium diboride bulk. Our findings for Au and Pb counter- 
electrodes were interpreted within the two-gap model for the energy spectrum of 
MgB2 and thus provide novel arguments in favor of a multi-band scenario for mag- 
nesium diboride. 



I 



I 



2, 2008 6:7 WSPC/INSTRUCTION FILE ingb2-MPLB 



Point- Contact Spectroscopy of Magnesium Diboride With Different Counter- Electrodes 11 
Acknowledgements 

This work was partially supported by the Polish Committee for Scientific Research 
(KBN) under Grant No. PBZ-KBN-013/T08/19 and by INTAS under project No. 
2001-0617. 



References 

1. C. Buzea and T. Yamashita, Supercond. Sci. Technol. 14, R115 (2001). 

2. H. Kotegawa, K. Ishida, Y. Kitaoka, T. Muranaka and J. Akimitsu, Phys. Rev. Lett. 
87, 127001 (2001). 

3. S. L. Bud'ko, G. Lapcrtot, C. Petrovic, C. E. Cunningham, N. Anderson and P. C. 
Canficld, Phys. Rev. Lett. 86, 1877 (2001). 

4. N. L. Bobrov, P. N. Chubov, Yu. G. Naidyuk, L. V. Tyutrina, I. K. Yanson, W. N. 
Kang, H.-J. Kim, E.-M. Choi, S. I. Lee, in New Trends in Superconductivity, ed. J. 
F. Annett and S. Kruchinin ( Kluwer Academic Publishers, Dodrecht, 2002), vol. 67 , 
p. 225. 

5. Yu. G. Naidyuk, I. K. Yanson, O. E. Kvitnitskaya, S. Lee, S. Tajima, cond-mat/0211134 
(2002). 

6. I. K. Yanson, V. V. Fisun, N. L. Bobrov, Yu. G. Naidyuk, W. N. Kang, E.-M. Choi, 
H.-,J. Kim, S.-I. Loo, cond-mat/0206170 (2002); I. K. Yanson, Yu. G. Naidyuk, O. E. 
Kvitnitskaya, V. V. Fisun, N. L. Bobrov, P. N. Chubov, V. V. Ryabovol, G. Behr, W. 
N. Kang, E.-M. Choi, H.-J. Kim, S.-I. Lee, T. Aizawa, S. Otani and S.-L. Drechsler, 
cond-mat/0212599 (2002). 

7. P. Samuely, P. Szabo, J. Kacmarcik, T. Klein, A. G. M. Jansen, cond-mat/0211544 
(2002). 

8. A. I. D'yachenko, V. Yu. Tarenkov, A. V. Abal'oshev, S. J. Lewandowski, cond- 
mat/0201200 (2002). 

9. J. M. An and W. E. Pickott, Phys. Rev. Lett. 86, 4366 (2001). 

10. Y. Kong, O. V. Dolgov, O. Jepsen and O. K. Andersen, Phys. Rev. B64, 020501(R) 
(2001). 

11. K.-P. Bohnen, R. Heid and B. Renker, Phys. Rev. Lett. 86, 5771 (2001). 

12. A. Y. Liu, I. I. Mazin and J. Kortus, Phys. Rev. Lett. 87, 087005 (2001). 

13. H. J. Choi, D. Roundy, H. Sun, M. L. Cohen and S. G. Louie, Phys. Rev. B66 
020513(R) (2002). 

14. H. J. Choi, D. Roundy H. Sun, M. L. Cohen and S. G. Louie, Nature 418, 758 (2002). 

15. H. Schmidt, J. F. Zasadzinski, K. E. Gray and D. G. Hinks, Physica C385, 221 (2003). 

16. J. Kortus, I. I. Mazin, K. D. Belashchenko, V. P. Antropov and L. L. Boyer, Phys. 
Rev. Lett. 86, 4656 (2001). 

17. S. V. Shulga, S.-L. Drechsler, H. Eschrig, H. Rosncr, W. E. Pickett, cond-mat/0103154 

(2001) . 

18. I. I. Mazin, O. K. Andersen, O. Jopscn, O. V. Dolgov, J. Kortus, A. A. Golubov, A. 
B. Kuz'menko, D. Van der Marcl, Phys. Rev. Lett. 89, 107002 (2002). 

19. A. Brinkman, A. A. Golubov, H. Rogalla, O. V. Dolgov, J. Kortus, Y. Kong, O. 
Jepsen, O. K. Andersen, Phys. Rev. B65, 180517(R) (2002). 

20. G. E. Blonder, M. Tinkham and T. M. Klapwijk, Phys. Rev. B25, 4515 (1982). 

21. H. Schmidt, J. F. Zasadzinski, K. E. Gray and D. Hinks, Phys. Rev. Lett. 88, 127002 

(2002) . 

22. E. L. Wolf, Principles of Electron Tunneling Spectroscopy ( Oxford University Press, 
New York, 1985). 



February 2, 2008 6:7 WSPC/INSTRUCTION FILE ingb2-MPLB 



12 A. I. D'yachenko et al. 



23. R. S. Gonnelli, A. Calzolari, D. Daghero, G. A. Ummarino, V. A. Stepanov, G. 
Giunchi, S. Cerasara and G. Ripamonti, Phys. Rev. Lett. 87, 097001 (2001). 

24. X. Z. Liao, A. Serquis, Y. T. Zhu, J. Y. Huang, L. Civale, D. E. Peterson, F. M. 
Mueller and H. F. Xu, cond-mat/0212571 (2002). 

25. O. N. Dorokhov, Pis'ma Zh. Eksp. Tear. Fiz. 36, 259 (1982) [JETP Lett. 36, 318 
(1982)]; O. N. Dorokhov, Solid State Commun. 51, 381 (1984). 

26. K. N. Schep and G. E. W. Bauer, Phys. Rev. Lett. 78, 3015 (1997); K. N. Schep and 
G. E. W. Bauer, Phys. Rev. B56, 15860 (1997). 

27. Y. Naveh, V. Patel, D. V. Averin, K. K. Likharev and J. E. Lukens, Phys. Rev. Lett. 
85, 5404 (2000). 

28. A. Bardas and D. V. Averin, Phys. Rev. B56, R8518 (1997). 

29. R. C. Dynes, V. Narayanarnurti and J. P. Garno, Phys. Rev. Lett. 41, 1509 (1978). 

30. J. G. Simmons, J. Appl. Phys. 34, 238 (1963). 

31. A. B. Kuz'menko, F. P. Mena, H. J. A. Molegraaf, D. Van der Marel, B. Gorshunov, 
M. Drcsscl, 1. I. Mazin, ,J. Kortus, O. V. Dolgov, T. Muranaka and J. Akimitsu, Solid 
State Commun. 121, 479 (2002). 



